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Abstract: Multiple interactions between parton showers and the surrounding QCD matter are 
expected to underlie the strong medium-modifications of jet observables in ultra-relativistic heavy 
ion collisions at RHIC and at the LHC Here, we note that such jet-medium interactions alter 
generically and characteristically the color correlations in the parton shower. We characterize 
these effects in a color-differential calculation of the medium-induced gluon radiation spectrum to 
first and second order in opacity. By interfacing simple branching histories of medium-modified 
color fiow with the Lund hadronization model, we analyze how the medium modification of color 
correlations can affect the distribution of hadronic fragments in jets. Importantly, we observe that 
jet-medium interactions give rise to the medium- induced color decoherence of gluons from the parton 
shower. Since hadronization respects color fiow and since each color singlet in a parton shower is 
hadronized separately, this medium-induced color decoherence leaves characteristic signatures in 
the jet fragmentation pattern. In particular, it can contribute to the quenching of leading hadron 
spectra. Moreover, it can increase strongly the yield of soft hadronic fragments from a jet, while the 
distribution of more energetic hadrons follows naturally the shape of a vacuum-like fragmentation 
pattern of lower total energy. 
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1 Introduction 

In ultra-relativistic nucleus-nucleus collisions at RHIC and at the LHC essentially all hadronic par- 
ticle distributions at high transverse momentum show strong modifications if compared to baselines 
established in proton-proton collisions. In particular, one observes a strong suppression of all single 
inclusive hadron spectra up to the highest transverse momenta {pt ~ 100 GeV) analyzed so far 
[1, 2], and significant modifications of jet structures up to jet energies exceeding 300 GeV [3-5]. 
Characteristic features of this jet quenching phenomenon include its dependence on centrality and 
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azimuthal orientation (both of which yielding information on the path-length of in-medium propa- 
gation), its dependence on kinematic variables including the center of mass energy of the collision 
and the transverse momentum of the hard process, its approximate independence on the produced 
hadron species (at least in the limited low-p^ range in which the latter has been measured to date), 
and the absence of quenching effect for high-p^ prompt photons and Z-bosons. The totality of 
these data from LHC and the data from RHIC [6, 7] motivates a dynamical picture of jet quenching 
according to which partons are produced in nucleus-nucleus collisions via high-momentum trans- 
fer processes at standard perturbative rates, but lose energy and branch differently in the dense 
QCD medium through which they propagate. Starting with the seminal works of Baier, Dokshitzer, 
Mueller, Peigne and Schiff (BDMPS) [8]and Zakharov (Z) [9] in the 1990s, a large number of par- 
ton energy loss calculations aim at formulating and exploring this jet quenching phenomenon in a 
QCD-based setup (for recent reviews, see [10-13]) 

In parton energy loss calculations, one generally considers a class of processes in which high- 
energy partons, produced in nucleus-nucleus collisions, interact via gluon exchanges with the sur- 
rounding QCD matter while branching. Color exchange between the partonic projectile and the 
QCD medium is intrinsic to such interactions. As a consequence, the color connection within a 
high-energy parton shower and between the shower and the rest of the event will be modified by 
the medium. Hadronization respects color correlations in mapping partonic color-singlet configu- 
rations into hadrons and can thus be sensitive to medium-modified color flow. Since the input of 
any hadronization routine is different if the color connections are different, medium modifications 
of the hadronization process may be expected to persist even if the latter occurs time delayed and 
thus outside the QCD medium. Taking such qualitative considerations into account, several groups 
have explored heuristic models of medium- modified hadronization in the recent past, pointing to 
possible changes in the hadrochemical composition of jets [14] or in the yield and distribution of 
baryons [15]. Also, some basic scenarios of medium-modified hadronization were explored in Monte 
Carlo simulations of medium-induced parton energy loss [16]. Notwithstanding these efforts, QCD- 
based calculations of jet quenching have remained mainly focused on kinematic changes of parton 
branching. Here, we extend the calculation of medium-induced gluon radiation of BDMPS-Z to 
the study of the color-differential case. We investigate, in particular, which medium-modified color 
connections can arise in the interactions of a parton shower with a QCD medium. Based on these 
calculations, we shall then test the response of hadronization models to a medium-modified color 
flow and discuss conceivable experimental signatures. 

Our work is organized as follows. In section 2, we set the stage by discussing qualitative features 
of the color flow of a parton shower developing in a QCD medium and how medium-modified color 
connections would affect the input of standard hadronization models. Sections 3 and 4 provide an 
explicit color-differential analysis of the medium-induced gluon radiation spectrum to, respectively, 
first and second order in opacity. We discuss which new medium-induced color correlations arise, 
how the weight of these contributions depends on the relation between the momentum and length 
scales in the problem (i.e. on the formation times), and how higher orders in opacity enhance the 
probability that the radiated gluons are decohered in color from the leading partonic fragments. 
The readers who wonder more about the conceivable effects of medium-modified color flow on 
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jet observables, but are less interested in the technical aspects of our calculation, may jump in a 
first reading directly to Section 5. There, we investigate how the output of standard hadronization 
models will change if they are interfaced with the characteristic medium-modified color flow patterns 
identified in our calculations of sections 3 and 4. Much of this discussion will focus on the Lund 
string fragmentation model, but we discuss also implications for cluster hadronization models. 
Sec. 5, read together with Sec. 2, provide a sufficiently self-contained digest of our qualitative 
arguments and of the essential message conveyed by our paper. We finally summarize and discuss 
our findings and provide an outlook to further open questions. 

2 Medium modification of parton splittings and their color connections 

To set the stage for the study of medium-modifications of high-p^^ processes to be carried out in this 
paper, we start from the 'vacuum' baseline of a hard partonic interaction in the absence of medium 
effects. To be specific, we first consider the case of a quark of color / from 'hadron 1' (proton/nucleus) 
that hard scatters on a quark of color i from 'hadron 2' (proton/nucleus) and subsequently radiates 
a gluon in the showering stage. To leading order in the coupling constant (when the hard 
scattering occurs via a single gluon exchange) and to leading order in the number of colors 
(when gluons can be represented as quark-anti-quark pairs), this vacuum process is depicted in 
Fig. 1. In order to ensure color neutrality, the colliding hadrons - here and in the following - 
will be depicted schematically as two opposite color charges. In the following sections, we shall 
calculate medium-modifications of such color-differential partonic processes and discuss how the 
medium-modified color flow affects their hadronization. Color connections are shown schematically 
in Fig. 1 by supplementing the quarks of color / and i in the incoming nuclei with anti-quarks of 
the corresponding color. These anti-quarks should be thought of as formal placeholders for what 
remains of a hadron once a quark of color I or i has been taken from it. They do not partake in 
the partonic interaction, but provide the color reservoir with which the scattered quarks are color 
correlated. 

The dynamics underlying hadronization is not understood from first principles. It is known, 
however, that the color flow of the underlying perturbative process is relevant and this is im- 
plemented in modern, phenomenologically successful models of hadronization, such as the Lund 
string-fragmentation in Pythia or the cluster-decay in Herwig. Here we discuss shortly how 
information about the color flow in the partonic process enters these models. 

Cluster hadronization models, as implemented e.g. in the Herwig event generator [17], group 
the result of a perturbative shower evolution into a set of color singlet clusters by splitting each 
gluon in the final state into a gg-pair; clusters are then decayed and hadronized independently. 
This is illustrated in Fig. 1 (left), where the quark k is combined with the anti-quark of the split 
gluon into a cluster. The most energetic hadron is then typically a fragment of this cluster, and 
the distribution of the fragments will depend on the momentum of the cluster and on its invariant 
mass. Alternatively, the Lund model, as implemented e.g. in the Pythia event generator [18], 
groups the same perturbative information into a set of color strings that start with a quark, follow 
the color flow by including gluons as kinks, and end on an anti-quark; these color singlet strings 
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Figure 1. The vacuum baseline: a hadron-hadron collision with a hard process followed hy a q ^ qg 
splitting in the shower stage. Color-flow in the large-A'^c limit is shown explicitly and the color-singlets 
(a cluster in the left panel, a string in the right panel) to be interfaced with the hadronization routine 
(Herwig cluster-decay or Pythia string-fragmentation, respectively) are displayed in color. 



are then hadronized (through excitation of qq pairs from the vacuum) according to a prescription 
that requires kinematic information about both the end-points of the string and all the kinks. The 
multiplicity and distribution of final hadrons will then depend significantly on the 'length of the 
string', that is on the separation of the quark and anti-quark end-points in momentum space. 

The presence of a medium with which high-energy partons can exchange color can clearly alter 
the color connections described above and, by changing the properties of the clusters/strings, have 
an effect on the final hadron spectra. In the remainder of this section we illustrate the essential 
ideas, focusing on the simple situation of a single interaction of the hard parton with the medium. 
We treat the cases of an incoming quark and gluon separately. 

2.1 Medium-induced color flow for a quark projectile to first order in opacity 

'Jet-quenching' calculations consider the interaction of high-energy partons in the dense QCD 
medium produced in heavy ion collisions. In Figs. 2 and 3, we display the simplest case of such an 
interaction: an elastic scattering of a high-energy quark in the plasma induces the radiation of a 
gluon. The interaction between the quark projectile and the medium occurs through the exchange 
of a single gluon, depicted here as a gg-pair in the large A'c-limit. 

For the configuration in Fig. 2, the 'leading' color singlet cluster, the one containing the quark 
k, shows the same color structure as the vacuum baseline, the color-flow connecting the high-pr 
quark with the anti-quark component of the gluon. In the LUND picture, the string containing the 
leading quark k links to the radiated gluon as in the vacuum baseline in Fig. 1, but now ends on 
an anti-quark from the medium. However, similarly to the vacuum case, this end point sits at low 
transverse momentum (with respect to the energy of the hard parton). For the above reasons, we 
shall refer to the color configuration depicted in Figure 2 as vacuum-like or - emphasizing the link 
of the radiated gluon with the projectile fragment - projectile- connected. 

Figure 2 is only one possible color structure that can emerge from a single interaction of the 
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Figure 2. Example of a hard qq qq event embedded in a nucleus- nucleus collision in which one of the 
high-p7^ quarks interacts once with the surrounding QCD matter which induces gluon radiation. Gluons 
are denoted by gg-pairs. The red lines denote the color singlet into which the leading quark k is grouped 
to form a cluster (left-hand side) or a LuND-string (right-hand side) in the corresponding hadronization 
models. 




Figure 3. Same as Fig. 2, but for the case that the high-pr quark radiates the gluon prior to interacting 
with the medium. 



projectile with the QCD medium. The second possibility is shown in Fig. 3, where (from the point of 
view of color flow) interaction with the medium occurs after the gluon emission. As a consequence, 
the leading color singlet cluster combines a quark at projectile energy with a target component 
at low (thermal) pr- In [19] the invariant mass of this cluster was shown to be parametrically 
larger than the one of the cluster in Fig. 2. Analogously, in the Lund framework the leading 
string connects the quark k directly to the target. The radiated gluon is, in both descriptions, 
color decohered from the projectile and will contribute only to an increase of the multiplicity of 
soft hadrons. In the following, we shall refer to these color configurations as medium-modified or 
gluo- decohered. 

We finally relate this discussion to the diagrams in Fig. 4 that are usually drawn for the 
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Figure 4. The {N = 1 in opacity) gluon-radiation amplitude arising from a single interaction of the 
projectile with the scattering centers of the medium. 



calculation of parton energy loss within an opacity expansion. In the figure, the initial hard process 
is sub-summed in a dark blob from which a single high-p^ quark line emerges at the light-cone time 
Xq (employing light-cone coordinates will be convenient in the following). The interaction of the 
high-pT- quark with the medium is described in terms of one-gluon exchange with a colored scattering 
center, or in terms of multiple such one-gluon exchanges. There is a direct correspondence between 
the Feynman diagrams in Fig. 4 and the ones displaying the color-flow in the large- iVc limit in Figs. 2 
and 3. One easily checks that the first amplitude in Fig. 4 corresponds to the color flow shown in 
Fig. 2, while the second diagram refers to the color configuration depicted in Fig. 3. Finally, the 
amplitude with the triple gluon vertex contributes to both color channels shown in Figs. 2 and 3 
(for details, see section 3). This correspondence prompts us to label the two color channels as Final 
State Radiation (FSR) and Initial State Radiation (ISR), depending on whether the gluon radiation 
occurs after or before the elastic scattering. 

A detailed color-differential calculation of medium-induced gluon radiation starting from the 
diagrams in Fig. 4 will be presented in Sec. 3. 

2.2 Medium-induced color flow for a gluon projectile to first order in opacity 

For the case of a high-py gluon produced in a hard scattering process, a larger number of color 
channels is involved. In Figure 5 we show the color configurations that arise to first order in 
opacity ^. For each case, we sketch only the associated Lund string that contains the most energetic 
gluon. 

In the large- iVc limit, a gluon is represented as a quark- ant i- quark pair. There are contributions 
to medium-induced parton branching in which both the scattering and the radiation occur off the 
q (q) leg of the projectile gluon while the q (g) leg is a silent spectator. In close analogy to the case 
of a quark projectile, we distinguish four cases in which - from the point of view of color flow - the 
gluon branching occurs either on the quark or on the anti-quark leg, and either after (Final State 
Radiation) or before (Initial State Radiation) color exchange with the medium. We refer to these 

^Hcrc and in the following, we distinguish the two daughter gluons in the medium-modified g g g splitting as 
leading (i.e. carrying 1 — ~ 1) and subleading (i.e. carrying Xg I), respectively. Exchanging these labels of the 
two daughter gluons would amount to exchanging the diagrams on the left and right column of Fig. 5. 
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four FSR(g), FSR(g), ISR(g), ISR(g). In addition, there are two channels in which the 

gluon radiation is connected, in the large- A'^c limit, to the q-leg (^-leg) of the high-p^^ gluon while 
color-exchange with the medium occurs via the g-leg (g-leg), see last row of Fig. 5. We anticipate 
that these last two channels will provide a subleading contribution to the spectrum in the case of 
soft {xg ^ 1) gluon radiation. An explicit evaluation of the color-differential radiation spectrum 
will be presented and discussed further in Sec. 3 and in appendix A. Here, we note that the Lund 
string including a high-p^ gluon will inevitably have both end-points at (anti) quarks that typically 
carry low transverse momentum. It is a particularly 'long' string in the sense that it stretches from 
low-pT to high-pT and back again. 



3 Color flow and radiation spectrum in the = 1 opacity expansion 

In this section, we provide a QCD-based color differential calculation of the medium induced gluon 
radiation process discussed in Sec. 2. We do so within the framework of parton energy loss first 
formulated by Baier, Dokshitzer, Mueller, Peigne and Schiff (BDMPS) [8]. The QCD medium is 
parametrized, following the Gyulassy-Wang model [20], as a collection of static scattering centers 
at a discrete set of space time points Xn giving rise to a potential A^{x). We work in light-cone 
coordinates and in a high-energy approximation in which (in the light-cone gauge A'^ = 0) only the 
'— ' light-cone component for the scattering potential is relevant 



A-(x) = 



N 
n=l 



dq 

(2^ 



U(q) 5{x' 



(3.1) 



Here, R denotes the representation of the parton suffering the elastic scattering. It is customary to 
formulate calculations of parton energy loss in a rotated frame in which the longitudinal axis points 
along the initial direction of propagation of the high transverse momentum parton. This parton 
has thus a large initial longitudinal momentum and no initial transverse momentum. During its in- 
medium propagation it will accumulate transverse momentum q from interactions with the medium, 
and will lose a light-cone energy fraction Xg by emitting a gluon of transverse momentum kg. For 
a parton of mass M and virtuality Q, the incoming and outgoing momenta read 



Pi 



Pf 



^ [q - kgf + ■ 
1 - Xg)p+, — q-kg 

2(1 - Xg)p + 



and for the gluon four-momentum and polarization vector one has 



(3.2) 



XgP^ 



' 2XgP~^ 



kn 



0, ■ 



XgP 



+ ' ^9 



(3.3) 



To illustrate the essential ideas within a sufficiently simple setup, we calculate in this section explicit 
color differential expressions for the gluon radiation at A^ = 1 order in the opacity expansion. 
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Figure 5. The six different color configurations tliat arise to first order in opacity for a fiigli-pr gluon. 
Curved colored lines represent the 'leading' Lund strings, i.e. the ones connected to the most energetic 
parton. Representing the high-pr gluon as a qq-pair, there are 4 contributions corresponding to scattering + 
Final State Radiation on the same quark or anti-quark leg (diagrams in first row) and to the corresponding 
terms for the Initial State Radiation case (diagrams in second row): in these four channels the second 
component of the gluon acts as a spectator. The contributions in the third row for which the radiated 
gluon is emitted from the q-leg while the medium couples to the ^-leg of the high-p^ gluon (or vice- versa) 
are suppressed in the soft {xg <^ 1) limit. 

3.1 Color-differential gluon radiation off a quark: the case = 1, x'^ = — oo 

We consider first the particularly simple case of a high-energy quark produced in the distant past 
{xq = — oo) and that interacts just once with the medium. The diagrams contributing to this 
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(a) (b) (c) 



Figure 6. The three diagrams contributing to the gluon radiation off an on-shell parton suffering a single 
scattering in an external color field. 



process are shown in Fig. 6. The amplitude corresponding to gluon emission after the scattering 
(first diagram in Fig. 6) takes the simple form 

N 

Pf^9 



n=l Pf''^9 



(3.4) 



Analogously, one can write the expressions for the other two Feynman diagrams shown in Fig. 6. 
The sum of these three contributions can be written in a compact form 



N 

-2ig 



n=l 



^a-{kg-q) 



k'g+xlM^ 



-X; 



!M2 



(3.5) 



Here and in the following, we work in the soft {xg <^ 1) limit in which terms sub leading in Xg are 
neglected. However, we have kept in Eq. (3.5) terms of order x^M^, so that the discussion can 
readily be extended to deal with the heavy-flavor case. The color-inclusive spectrum of radiated 
gluons is obtained from |A^p after averaging (summing) over the initial (flnal) states: 



1 , da'^'^ 

— rkZ 



CA-,{[K,-K,f) 



IT 



where the notational shorthands 



kr 



kj + xlM^ 



kr 



-xjM^ 



(3.6) 



(3.7) 



have been used, and the momentum kicks q received from the medium were averaged according to 
the corresponding elastic cross section: 

1 rf(7^'^ 



(...) 



dq 



a"^ dq 

In the M — > limit, the well-known Gunion-Bertsch spectrum 



dq\A{q)\' 



dNr, 



3 dkgdk+ 



C, 



a. 



kq{kg 



(3.8) 



(3.9) 
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is recovered. Notice that once the radiation spectrum is normahzed by its elastic cross section, the 
overall color factor Ca is universal, irrespective of whether the incoming projectile is a quark or a 
gluon. The elastic scattering cross section entering Eq. (3.8) is often chosen to be of Yukawa-type, 
|^(g)P oc ji\)/{q^ + /uf))^, but we shall not rely on a specific functional shape of |^(q)p in the 
following. 

Extending the derivation of the Gunion-Bertsch spectrum to the color-differential case is straight- 
forward. After identifying in the three-gluon vertex oc the two different 'color-orderings', 
one writes the total radiation amplitude as: 

iM = iM''''' ^iW^'' , (3.10) 

where A^""^ ocf^t"-^ and A^"i" cx:t"H", respectively. One then notes that in the cross section the 
interference between the two contributions is suppressed by a factor since TT{t°'t°'^t°'^t°-) = 

ClNc and Tr(rrirt"i) = -{1/2N^)CfNc. Thus, 



'''"'-^t + ^t + Oil/Nl), (3.11) 



from which one obtains, to leading order in Nc-, the spectrum of radiated gluons for each color 
channel. Keeping track of the exact dependence on the energy fraction carried by the emitted gluon 
one finds 



dNg 



^ die gdk^ 



(3.12) 



^ dkgdk + 



aai 



aia 



2 TT 

2 7r2 



Ko-K^Y), (3.13) 



where we have used the shorthand Kq = (fc ^x^qp^iM'^ ■ leading order in Xg (i.e. in the soft 
limit), one sees that both color channels have exactly the same weight and that the sum of the two 
contributions yields Eq. (3.6). In general, however, there is no general argument for why different 
color channels should contribute equally. To illustrate this point already here, we have kept in Kq 
the subleading dependence on Xg. 

3.2 Color-differential gluon radiation off a gluon: the case = 1, Xq = — oo 

So far, the calculations in this section have focused on the case of a projectile quark that radiates 
a gluon in response to medium-induced scattering. For an anti-quark projectile, the result is com- 
pletely analogous and does not require separate discussion. For a gluon projectile, however, one 
finds a larger set of distinct color configurations, as discussed already in section 2.2 and shown in 
Figure 5. Here we discuss this case explicitly for a gluon produced in the distant past. We first 
note that, identifying in Fig. 6 the projectile line with a gluon, three Feynman diagrams (a), (b) 
and (c) contribute to the radiation amplitude. To be specific, we label the color of the incoming 
gluon as d. This gluon splits into the most energetic 'projectile'-fragment of color b and a radiated 
gluon of color a. The color exchanged with the medium is labeled as ai. In terms of the generators 
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Ta of the adjoint representation, the corresponding contributions to the radiation amphtudes have 
the color structure: 



M(a) oc {TlTl%, , oc {n'n)M ^ -^(e) « {1^,^%^ ' (3-14) 

Expressing the generators Ta in terms of the structure constant of the Lie algebra, one can conve- 
niently exploit relations like {T"'T°'^)m = 2 Tr ([t'', t"] [t'^^, t'^]) to identify the color channels involved. 
Due to the invariance of the trace for cyclic permutations only 6 distinct color configurations arise. 
There is a one-to-one correspondence between these 6 independent permutations and the six color 
differential diagrams for medium-induced gluon radiation shown in Figure 5. One finds 

ISR(g) O haiad , ISR(g) ^ hdaai , 

FSR(g) o haaid , FSR(g) ^ hdaia , 

RgSg ^ hadai , RgSg ^ haida . (3.15) 

In the calculation of the radiation cross section, interference terms between these six contributions 
correspond to non-planar diagrams and are subleading in the large-A^c limit. In appendix A, we 
provide explicit expressions for the amplitudes (3.14) and their six distinct color-differential contri- 
butions (3.15). From these expressions, one sees easily that the contributions RgSg and RgSg are 
suppressed by a factor 0{xg) compared to the other four terms. Up to 0{xg) corrections (hence 
consistently with the approximation assumed in getting the color-inclusive GB spectrum), the terms 
that describe initial or final state radiation off a q or g-leg provide color-differential radiation cross 
sections of equal size. After summing over the final and averaging over the initial polarization of 
the gluons, each of these four contributions takes the form 



dk^dk 



" " ' ^ ^ (3.16) 



4 7r2 \ k\k-qy 



with 2 = FSR/ISR(g/g). Summed together, these 4 leading contributions lead to the inclusive result: 

,+_dN ^ a, / \ 

k zjznz = Ca— { — — ) ■ (3.17) 



dk+dk 7r^\k^{k-q 

In summary, these results show that at small Xg and in the large Nc limit, the = 1 radiation 
spectrum of a projectile gluon can be viewed as an incoherent superposition of the medium-induced 
radiation of a quark (with a silent anti-quark spectator), and of an anti-quark (with a silent quark 
as spectator). Modifications to this picture, arising from the qualitatively novel color configurations 
(RgSg and RqSq), are suppressed in the soft limit by a factor O(x^). 

3.3 Color-differential medium-induced gluon radiation: the case N = 1, Xq = 

The gluon spectrum (3.9) vanishes in the absence of a medium, consistently with the idea that a 
parton coming from the distant past is on-shell. The situation realized in a hadronic collision is 
different. There, high-energy partons are produced around the time of the collision, at x^J" = 
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Mo 

Figure 7. The vacuum-radiation diagram: tlie dark blob denotes the hard process in which the high- 
momentum parton is produced (off-sheU) inside the medium. In the right panel we keep track of the color 
flow. 



say, and they branch even in the absence of further in-medium interactions. Here, we calculate 
the medium-modification of this vacuum radiation. That is, we consider the N = 1 corrections 
in opacity to the 'vacuum' splitting process shown in Figure 7. The amplitude for this vacuum 
baseline is written as 

^Mo = (.^t")[(2p; + A:,)-6,] ^^^^^^^^^_^^ J(p^ + fc,)e-(^/+^-)-° , (3.18) 
where the current J{pf + kg) represents the hard production process. The resulting spectrum reads 

with its characteristic coUinear divergence in the massless limit and dead-cone effect suppress- 
ing the radiation of small-angle gluons in the case of emission off a massive quark, where = 
feg/ (fc^ -|- x^M^) . In general, the radiation amplitude can be expanded in the number of gluon 
exchanges with the medium, 

^A^rad = iMo + iMi + iM2 + ... (3.20) 

To perform a calculation that is complete to first order in opacity, the above expansion must include 
not only the complete set of one-gluon exchanges between partonic projectile and medium Aii, but 
also a specific subset, which we label A^™'^*, of two-gluon exchanges. Contributions to first order 
in opacity arise both from (|A^ip) and Re (A^2-A^o)- 

The medium average (...) involves a color trace over the target and an integration over the 
transverse position of the scattering centers. For the two scattering centers that enter the term 
Re(A^2-Mo); these averages reduce to 

Tr (t^:)T^^,^) = Snn'S'^^'^'TF/A, J dx^c-'^'^^^'^-^^- = 2'K5{q^ + q^) , (3.21) 

where Tp = (1/2) and = A^^ for scattering off a quark or a gluon from the medium, respectively. 
As a consequence, after average over the target, the non- vanishing contributions to Re {M.2J^o) arise 
from processes where the two gluons link to the same scattering center in the medium and where 
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I I I 

Figure 8. The gluon-radiation amplitude arising from interacting twice with the same scattering center of 
the medium. In the = 1 opacity expansion, after the target average (over color and Xn), the diagrams 
contribute to the single-gluon spectrum with an interference term decreasing the weight of the vacuum 
radiation. 



neither color nor transverse momentum is transferred from medium to the fragmenting q q + g 
system. It is usual to refer to these contributions as contact terms and they are needed for a 
formulation that conserves probability, thus playing the role of virtual corrections. The contact 
terms contributing to A^2 are shown in Figure 8. To first order in opacity, one has 

(|Alo + 7Wi+7Wr* + -l'> = |-Mo|' + (|-Mi|'> + 2Re(A^r*A^*) + .... (3.22) 

The radiation amplitude can be recast in a form more suitable for a color differential calculation: 

M0 + M1+ = -M''"^ + M''''' + A^" , (3.23) 

where 

Ml = M"""' + M''''' and Mo + = M". (3.24) 

In contrast to the left-hand side of (3.23), interference terms between different contributions on the 
right hand side vanish to leading order in Nc. This allows us to derive, in the large-iVc limit, the 
color-differential radiation cross sections by squaring the individual contributions from the various 
color channels. More explicitly, the M"'"'^ term reads 

n=l ^ ' 

Ko-[l- e-^^i(^*-^o+) j Ki . (3.25) 
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In this expression, the term proportional to Kq corresponds to a hard parton that emits a gluon 
after scattering; it is the term A4i(a) in Fig. 4. On the other hand, the three-gluon vertex in the real 
emission amplitude A^i(c) in Fig. 4 gives rise to a term proportional to and the f^f^^-part 

of this commutator provides the contribution oc Ki in Eq. (3.25). The two Ki terms (with their 
different phase factors) arise from processes in which the propagator of the hard projectile parton 
is on or off-shell, respectively. In close analogy, one finds 

where the term proportional to Kq corresponds to A^i(fe) in Fig. 4 and the term proportional to 
Ki is the part of Aii(c) proportional to t°'^t°'. Again, the phase factors are different depending on 
whether the projectile from the hard event is on or off-shell. In the color-differential amplitude 
(3.26), the phase factors are written in terms of the 'transverse energies' 

Uq = — , Ui = . [6.Z( 

In general, the inverse transverse energies I/cJq, 1/uJi act as formation times for medium-induced 
gluon emission. The color-averaged result for medium-induced gluon radiation to first order in 
opacity is known to depend only on l/cJi (we reproduce this result as a check in Eq. (3.31) below). 
For the color-differential case studied here, we observe that the result depends on both formation 
times 1/cJo and 1/cJi. In particular, after squaring the above amplitudes and averaging over the 
longitudinal position x^ G [xq] Xq + L+] for a medium of constant density, we find a vacuum-like 
contribution from the color-channel 'aai' and a medium-modified one from the channel 'aia': 

(l-Mrr) - ((^0 - + Ki + 2K,.{K, - i^i)^^^^) , (3.28) 

(m-^2(:-?^^)KS.2((:-^jK; 

_ smKL^I _ s,n|.,.L^| ^ s,u|(^. -c.o)L^| W„ . jf. V (3.29) 

V ujqL+ uJiL+ {uJi-uJo)L+ J I 

To first order in opacity, the calculation is completed by including the channel 'a' in which no color 
is exchanged with the medium. We have (|A^^p) = |A^oP + 2Re (A^2°°*-^o)' where |A^oP accounts 
for the vacuum branching and 



2Re {M,Ml) ~ -K^ - 2 1 - ^^^^^ Kl 



sin[a;oi^^ 
sin[cZ;o-^'^] sin[(aJi — 'Ujq)L^\ 



Ko-KA . (3.30) 
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The sum (lA^^^'p) + (|A^f p) + 2Re (A^2°°*-^o) combines to the color- averaged medium-induced 
gluon radiation cross section 



dl 



med 



dk~^ dkn 



■el 



-"■0 



Kl)[l- 



(3.31) 



This expression shows exphcitly that the opacity expansion is an expansion in powers of L^/A^, 
where denotes the elastic mean free path of a gluon and R is the representation of the projectile. 
Here we are considering a quark, hence Cr = Cf- 

3.4 Dependence of color-differential = 1 gluon radiation on formation time(s) 

In the color-averaged cross section (3.31), the phase factor 

_ sin {u^ L+) \ _ r for l/wf > L+ , 
L+ j ~ I 1 for l/wf < L+ . 



(3.32) 



interpolates between the coherent and incoherent regimes. Thus, the color-inclusive medium- 
induced spectrum in Eq. (3.31) vanishes if the gluon formation time prior to its rescattering l/cuf is 
larger than the in-medium path length. This result matches the naive expectation that the medium 
can modify only the radiation of those gluons that are fully formed within its finite extension. 
Here, we discuss the different behavior of the various color channels: vacuum-like (aai), medium- 
modified (flia) and not color-correlated with the medium (a). They sum up to the inclusive result 
dl = dr^"^ + dl'^^'^, where 



k' 



dp 



dk^dkr, 



k' 



dP 



dk+dkn 



and 



k' 



dl 



med 



dk^dkr. 



k' 



dl 



med 



+ k+ 



dl 



mod 



+ k^ 



dl 



med 



aia 



dk~^ dk f] 



(3.33) 



(3.34) 



aai 9 

For the color-differential contributions on the right hand side of this equation, the in-medium path 
length L"*" has to be compared with the two distinct formation times I/Uq and 1/cJi. To discuss 
this dependence, we focus on the four limiting cases in which either one or both formation times are 
either significantly larger or significantly smaller than L"*". In order to display more transparently 
the physical meaning of the various terms we will exploit the large- Ai'c identities Ca = '2Cf for the 



color factors and A|j" 



2Aj" for the elastic mean free paths. 



1. Totally incoherent case {L~^ ^ 1/wj, with i = 0,1) 

In the limit — )■ oo, the color-averaged medium-induced gluon radiation spectrum (3.31) is 
proportional to {{Kq — Ki)^ — -\- K^y This is the incoherent superposition of: the usual 
Gunion-Bertsch spectrum for gluon emission oc {{Kq — Ki)'^y, a negative contribution to the 
vacuum radiation spectrum, oc — {K"^, that corrects the vacuum term for the probability 
that the radiated gluon interacts with the medium; and a vacuum-radiation spectrum that is 
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shifted in transverse momentum due to the rescattering of the radiated gluon in the medium 
oc (-K'l)- For the corresponding color differential contributions, we find from (3.28), (3.29) 
and (3.30) 



dl 



mod 



dl\j~^ dk f] 



dl 



mod 



dk'^dkn 



aia 



{{{K,-K,f) + {K\)) 



L+ 



CpKl 



k- 



dl 



mod 



dk^dk„ 



(3.35a) 
(3.35b) 
(3.35c) 



After inspection of the color factors and of the mean-free-paths involved, the above terms 
admit a transparent physical interpretation. Eq. (3.35a) is the sum of half of the GB spectrum 
{{Kq — KiY) by an on-shell quark and half of the contribution of reshuffled vacuum radiation 
{K\) by an off-shell quark. Eq. (3.35b), on top of this, gets an additional contribution from 
the vacuum radiation i^Q by an off-shell quark which then suffers a further elastic scattering 
in the medium. In this limit, the medium-modified color differential contribution 'aia' is the 
largest. As a consequence, the radiated gluon is decorrelated in color from the leading partonic 
fragment in more than half of the medium-modified parton branchings. This case is depicted 
for instance in Fig. 3. 

Totally coherent case (1/cJj ^ L^) 

For very large gluon formation times, the gluon is produced far outside the medium. At the 
level of the color-inclusive result the medium is not a source of an enhanced rate of parton 
branching: the color-averaged medium-induced spectrum in Eq. (3.31) vanishes. However, 
with the probability that an elastic interaction occurs, the fragmenting projectile is color 
connected to the medium ('aai' channel) rather than to the hard production process ('a' 
channel) : 



k' 



dl 



mod 



dk^dkn 



aai 



dl 



mod 



dk'^dkn 



0. 



aia 



k- 



dl 



med 



dk'^dkr, 



w,L+^0 A+ Tl^ 



(3.36a) 
(3.36b) 
(3.36c) 



The rearrangement between the 'aai' and 'a' channels reflects the fact that color is exchanged 
between medium and projectile even if the color averaged spectrum (3.31) remains unchanged, 
i.e. even if there is no overall enhanced probability of gluon radiation induced by the medium. 
Notice, however, that in both channels the gluon is color correlated with the highest-py 
fragment. 
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3. 1/cJi > L+ > l/cJo 

In this kinematic range, the color-inclusive medium- induced radiation vanishes. Remarkably, 
however, the color connection amongst the most energetic fragments is still medium-modified 
and two-thirds of the gluons radiated due to medium effects (and compensated by a corre- 
sponding depletion in the vacuum spectrum) are carried by the 'aia' contribution and are 
thus decorrelated in color from the most energetic fragment 

L+ 



dk~^dkg 



—Cp^Kl (3.37a) 



dl 



med 



dk~^dkg 



k' 



dl 



med 



dk~^dkg 



2^CA^KI , (3.37b) 
^o^:j-°^_s£^Cf^KI. (3.37c) 



The factor 2 difference between the 'aai' and 'aia' channels reflects the fact that Eq. (3.37a) 
gets contribution only from processes in which the propagator of the hard quark is on-shell 
prior to scattering and radiation; for contributions to Eq. (3.37b) this propagator can be either 
on- or off-shell. 

4. 1/ZUo > L+ > 1/Ui 

Also in this case, the medium-modifled color channel 'aia' - in which the gluon decorrelates 
in color from the leading parton - tends to be the most likely one for parton splitting 

~ Tr^ri\(^o"^i) +-^^1/ ' (3.38a) 



k' 



dk~^dkg 
, dl'^^'^ 

k 



dk~^dkg 

(i/med 

k^ 



dk~^dkg 



TUiL+^oo A+ 2 71- 
aai g 

(3.38b) 

aia g 

":r"°^^5(-^o-2Ko.Ki>. (3.38c) 

a ^ g 



As an aside, we note that there are kinematic arguments for why this limiting case may be 
less relevant for a medium-modifled parton shower^. 

In summary, whenever medium-induced gluon radiation is sizable, the qualitatively novel medium- 
modifled color connection 'aia' of the projectile with the medium arises in more than half of the 
processes. Moreover, even in cases in which the color-inclusive induced spectrum (3.31) is negligible, 
the presence of the medium can act as a source of color-decorrelation between the radiated gluon 
and the leading parton. 



^We recall that the gluon in the final state has transverse momentum kg and therefore a formation time l/oJo- 
The term 1 /wi can thus be viewed as the formation time of a gluon that did not yet undergo scattering with the 
medium. Since gluon emission is dominated by collinear branching, the initial gluon transverse momentum kg — q 
will be small in most branching processes, and a further momentum transfer from the medium is more likely to 
increase the gluon's transverse momentum than to reduce it. So, on qualitative grounds, the ordering \kg — q\ < \kg\ 
(and equivalently I/cjq < l/<^i) is more likely to occur than the opposite one. 
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4 Higher orders in opacity: the N = 2 case 

For a parton shower that develops in the vacuum, the two daughters of a. q ^ qg or g ^ gg 
branching process are always color connected. In contrast, in the = 1 opacity expansion carried 
out in the previous section, one finds that the interaction with the medium during a branching 
process decoheres the daughters in exactly half of the cases for a projectile parton produced in the 
distant past {xq = — oo), and in more than half of all cases whenever there is a medium-induced 
gluon radiation from a parton produced in the medium [xq = and 1/uJi ^ L~^)- Here, we discuss 
how these findings change once more than one gluon exchange between the partonic projectile and 
the medium is taken into account. 

In color-inclusive calculations, parton energy loss is largely determined by the average squared 
transverse momentum per unit path-length irrespectively of whether the momentum is transferred 
in one or several gluon exchanges. The fraction of radiated gluons that are color decorrelated from 
the most energetic projectile fragment may, however, be expected to increase with the number of 
interactions with the medium. Here, we support this expectation by performing a color-differential 
calculation of the N = 2 opacity contribution to gluon radiation by a parton produced in the distant 
past. 

As usual, the starting point is the expansion of the radiation amplitude in the number of gluons 
exchanged with the medium: 

M'""^ = Mo + Mi + Mf' + * + Mf' + MY' + ... . (4.1) 

As was already the case at A^ = 1, the above amplitude involves 'virtual' contact terms in which two 
gluons are exchanged between the projectile and a single scattering center in the medium, without 
net exchange neither of color nor of transverse momentum. Up to second order in opacity - when 
the projectile arrives on-shell from the far past, so that Aio = - the radiation spectrum receives 
contributions from 

|X™d|2 _ |_y^^|2 ^ |_y^dir|2 ^ \Ml"'\^ + 2Re Mf' M*, + 0{L/\,xf . (4.2) 

Here, a contribution to zeroth order in opacity is absent for a projectile coming from the distant 
past, since a radiated gluon cannot be emitted on-shell without interaction with the medium. For 
the same reason, interference terms 2ReA^2"^*-^o vanish in |A^'''^'^p. However, the contribution 
|A^2"^*P does not vanish since a color-neutral two-gluon exchange between the projectile and the 
medium can - in general - still transfer longitudinal momentum, allowing the radiation of a gluon. 
Thus the amplitude for a color-neutral two-gluon exchange with no transfer of momentum q with 
one scattering center, paired with a corresponding term in the complex conjugate amplitude, leads 
to a finite contribution. 

4.1 Direct contributions to A^ = 2 

To second order in opacity, the contributions involving two finite momentum transfers are depicted 
in Figure 9. Here and in what follows, we represent effective non-local gluon emission vertices as dark 
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Figure 9. The three gluon-radiation ampHtudes involved in the N = 2 scattering case. 



blobs. They denote the combination of the three radiation amphtudes in Fig. 6 and correspond 
to the Lipatov vertices in BFKL calculations. Their contribution to the radiation of a gluon of 
momentum kg after momentum transfer q from the medium can be expressed in terms of the 
currents ^ ^ 

^ k] + l^M^ ~ [kg - qf + x^M^ • ^^-^^ 
We also introduce the notational shortcuts 

Ji = J{kg,qi), J2 = J{kg,q2) and J12 = J{kg - q2, qi)- (4.4) 

Employing the above currents it is possible to express the N = 2 direct contribution in the conve- 
nient form suggested by Fig. 9 

Mf' = Mm + M[2]+Miu], (4.5) 

with 

iM[i] ~^t"2[r,t''^]2eg-Jie*^«^^ (4.6a) 

iM[2] ~^[t",t"^]r^2eg-J2e*^«^2 (4.6b) 

iMii2] ~^[[t^t'^^],^"^]2e3■Jl2e^^2^^e^(^o-^2)^^. (4.6c) 

It follows from the color structure of the above amplitudes that the N = 2 direct contribution 
can be organized into five distinct color channels 

[a2aai] , [a2aia] , [00201] , [01OO2] , [oi 02 o] . (4.7) 

An example is shown in Fig. 10. The vacuum-like contribution, in which the color flows in the 
large A^c-limit from the projectile quark to the gluon, corresponds to the term [002O1]. In the 
four other color channels, color flows from the leading quark projectile directly to the medium 
without passing through the gluon, i.e., the gluon is color decohered from the quark. To leading 
order in Nc, these five color channels do not interfere. From Eqs. (4.6a-4.6c) one readily obtains 
the radiation spectrum in the different color channel (interference terms being suppressed by 
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[ajaai] I > *i 




Figure 10. An example of a Feynman diagram contributing to gluon radiation at = 2 order in opacity, 
with the corresponding color-flows associated to it. Two different processes are accounted for: the radiation 
of a gluon from the second scattering center, after the incoming on-shell quark already suffered an elastic 
collision ('aa2ai' contribution to A^[2]) and the radiation of a on-shell gluon from the first scattering center, 
which then suffers a further elastic collision ('02001' contribution to A^[i2])- 



factors). Transverse momentum kicks from the medium are weighted by the corresponding elastic 
cross-section, employing the shorthand notation 



{{■■■))qi,Q2 — 




Furthermore, we introduce the current |J_ooP = {'^P^)'^ describing the parton coming from the far 
past. The vacuum-like contribution takes the form (after averaging over the longitudinal location 
of the scattering centers) 

i\MT'"f} - (1)^ |._.|^ {{iM' ^ ilJ.I^ . 2l^0flj..J..))^^^^ (4.9) 

Amongst the four other color channels, three do not carry phase interference terms, 

= (lAl^p) = ig^Cp (^^^ |J_ooP ^<|^|Ji2p^^ , (4.10) 
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while the fourth one has a shghtly more comphcated structure 



X+ J ' \\2' 2' 2' (woL- 



-2 ' J.. J.. +2 ' > J.. J..)) , (4,11) 

The sum of the above five terms yields the N = 2 direct contribution to the spectrum: 



|2) =Ag'CF (^] I J.ool" ( ( I Jil' + \J2\'' + 2| Ji2r - 2 ^ ,r7"V JiJ 



17 12 //i T 12 , I T 12 , 01 T 12 ~ COS(WO^"^) 



+ ^2 



2 



-2 ' - i°^';;'f ' J. . J. + 4 ^ - r,?-^"' J. \ . (4.12) 

(^02^+)2 (W2^ + )2 //<7i,q. 

From this and the corresponding color-differential expressions, one then checks easily that, in the 
totally incoherent limit, the vacuum-like contribution 

{\MT''''\'') _ ~ ^(|J2|'+|Ji2|'> (4.13) 

is always more than a factor 2 smaller than the sum over all five color channels 

(|Mf H _ ~ (I Jip + I + 2| Ji2p> . (4.14) 

This is consistent with the naive expectation that a vacuum-like fragmentation pattern is less likely 
to survive in higher orders in opacity. 

4.2 Virtual contributions to = 2 

The N = 2 opacity calculation is completed by the computation of the relevant contact terms. In the 
case at hand, a projectile arriving on-shell from the far past, these correspond to the terms (| A^2^''*p) 
and 2Re(A^3^^*A^2) Eq. (4.2). From the color point of view they contribute, respectively, to the 
'a' and 'aai'/'aia' channels: 

(l-Mff) = (l-Mff)a, 

2Re{Mt''Ml) = 2Re{Mt''Ml)aa, + 2Re{Mt''Ml)a,a. (4.15) 

The term (|A^2"^*P)a accounts for processes in which neither color nor transverse momentum is 
exchanged between projectile and medium, but where the longitudinal ('-' in light-cone coordinates) 
momentum transferred by the medium is sufficient to open the possibility of gluon radiation for a 
quark initially on-shell. This is a novel contribution, absent at N=l opacity. On the other hand, 
the terms 2Iie{Ai^"^Ail)aai/aia are virtual corrections which combine with the direct amplitudes 
and (|A^'f^°p), and that decrease the weight of the corresponding color channels. 
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The calculation of {\Ai™^\'^)a requires the evaluation of the same contact terms, shown in Fig. 8, 
as in the = 1 case for a quark produced at Xq = 0. Here, however, only the terms referring to a 
quark arriving on-shell are retained 

In the color inclusive spectrum, this contribution cancels the cUo-dependent terms in (4.12). We 
further note that most 'virtual' corrections in the opacity expansion can be regarded as proba- 
bility conserving terms that subtract yield from terms of lower order in opacity. In contrast, the 
contribution (4.16) is always positive, and it opens a channel 'a' absent at lower order in opacity. 




Pi 



CD 




Pf 



Figure 11. An example of the contact diagrams contributing to A^™*. 



The calculation of the term 2Re(A^3'''*A^*) involves the evaluation of diagrams like the one in 
Fig. 11. Details of the calculation are given in Appendix B. One gets for the two channels 



|. ,2// 3,^,2 1, ^ ,2 , ^1 -cos(a;o^+) 



'2 



A+y ' ' \\ 2' 2' {uJoL 



•31. '72 



and 



3|T|2 liTi2 cos(ci;o^v"'") ^ ^ , o 1 ~ cos{uj2L^ 



^iJip-^iJ2r-2 ^ ,rrZ:; ' JvJ2+2 ^ ;zri.: ' J2-Ji2)) . 



2' ^' (cJoL+y ' ' {cJ2L+y 

4.3 Gluon radiation up to order A^ = 2 in opacity 

We can now write the medium-induced radiation spectrum up to order = 2 in opacity. It can be 
written as dl = dl^^^ + dl'^^'^, where dl'^^^ is the lowest-order (A^ = 1 opacity) result: 
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while rfJ^=2 arises from the sum {\Mf\'^) + {\MY^\'^) +2Re{Mf'^Ml) . In the large- A^'c limit when 
Ca = 2Cf for the color charges and Ag = 2\g for the mean free paths, one finds 

N=2 - / r-^\ 2 



k- 



dl 



dk^dk 



Ca 



a. 



'12 



Jir + 2 — — 



1 — cos(a;2i^''') 



(Cl^2^+)^ 



2J2 J 



12 



(4.20) 



Qi,Q2 



which coincides with the well known result quoted for instance in Ref. [16]^. 

While the above color-inclusive result depends on a single formation time I/CJ2, the color- 
differential contributions to second order in opacity depend already on three different formation 
times and have a somewhat more complicated structure. Let us now express the radiation spectrum 
in a color-differential way. 

The channels 'aai' and 'aia', already present at order = 1, receive in the N = 2 calculation 
a correction from the virtual terms 



dl. 



N=2 



aai 



dr'=' + dL„ 

aai aai 



In the above dl^^^ ^ 



dl^=^ 



dl, 
{l/2)dl 



N=2 



aia 
N=l 



diZT + diaia 



(4.21a) 
(4.21b) 



while the N = 2 corrections can be obtained from 



Eqs. (4.17) and (4.18) and read: 



dl 



N=2 



dk+dk 



aai 



, ^ 1 -cos(u;o2i:+) 

+^ 7= — TTT^ J 1-^12 



1 — cos{ujoL~^) 
1 — cos(aJ2-L 



(P2L+] 



J1J2 

J 

^2 «^12 



(4.22) 



91:92 



dl 



N=2 



dk~^dk 



aia 



^ 1 — cos(a;oi^^ 



-JIJ2 



p i -C0S(C^2^+) J J 



. (4.23) 



<5ri.'J2 



Furthermore the N = 2 calculation open new channels absent at lower order in opacity. The 
spectrum in the 'a'-channel, involving no color exchange, follows from Eq. (4.16): 



k- 



dl 



N=2 



dk^dk 



(4.24) 



91,92 



Concerning the channels arising from Ai"^^^ one has, for instance, the channel with vacuum-like color 
connection of the leading fragment with the radiated gluon. 



k- 



dl 



N=2 



dk+dk 



a 



aa2ai 



— \J2\ + 7:l'^12l +2 



1 — C0s{ljJ2L~ 

ip2L+Y 



-J2'Jl2 



(4.25) 



9i,92 



^Note that in [16] Minkowski coordinates are employed and the mean-free-path Xq is replaced by 2Xg 
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The spectrum in the other channels can be easily obtained from the corresponding squared ampli- 
tudes quoted in Sec. 4.1. 

Here, we end by considering the totally incoherent limit in which is much larger than all 
formation times in the problem. This is the dominant phase space region for gluon radiation in 
the sense that gluon radiation is not suppressed by destructive quantum interference. One checks 
easily that, in this totally incoherent limit, the totally virtual term (|A^2^''*p) (the one contributing 
to the 'a'-channel) vanishes. The vacuum-like 'aoi' and medium-modified 'oia' probability con- 
serving factors 2Re(A^3'''*A^^) have equal negative weight oc — 3| Jip/2 — | J2P/2. However, for the 
terms with two finite transverse momentum transfers from the medium, (|Al2''^P), the contributions 
with medium-modified color-flow (the ones in which the gluon is color-decohered) have a weight 
oc (2| Jip -|- I J2P + 3| Ji2n that is larger than the vacuum-like one oc (| + 1-^121^)- This illustrates 
with an explicit calculation the general expectation that the weight of vacuum-like contributions to 
the total radiation spectrum decreases as higher orders in opacity are accounted for. 

4.4 General features of a color-differential analysis of medium-induced gluon radiation 

We now summarize what we learned from the technical analysis performed in the previous sections, 
in which we calculated medium- induced gluon radiation in a spatially extended QCD medium 
within an opacity expansion. The color-inclusive result is consistent with previously known expres- 
sions. In addition, our calculation provides for the first time color-differential information about 
medium-induced gluon radiation. We distinguish in particular vacuum-like contributions (in which 
the leading final state parton is color correlated with the radiated gluon) from medium-modified 
contributions (in which the radiated gluon is decohered in color from the projectile) ^. Up to 2nd 
order in opacity, we found three vacuum-like contributions 

[a] , [aai] , [00201] . (4.26) 

Here, the notation [a...] indicates contributions to the gluon radiation amplitude in which the 
SU{Nc) generators carrying open indices Cj appear in the order specified in the brackets. In the 
large- iVc limit, interference terms between different contributions to the amplitude vanish, and the 
notation [a ...] provides an efficient labeling of the contributions to the gluon radiation cross section. 

The configuration [a] arises for processes in which the final q-g pair is in the same (total) color 
configuration as the incoming quark. For a quark produced inside the medium, this includes the 
entire 0-th order opacity contribution, as well as a negative probability conserving term to order 
= 1. For a quark produced in the distant past, energy-momentum conservation ensures that the 
contribution [a] is absent to orders = and = 1 in opacity, but it is found to contribute to 
order N = 2 with positive weight. 

The other vacuum-like color configurations [a ai] , [a a2 ai], and more generally configurations of 
the form [a an a„_i . . . ai], arise first to n-th order in opacity as real terms with n non- vanishing mo- 
mentum and color exchanges between medium and projectile parton. Higher orders in opacity will 

^More precisely, we call the large- -/Vc limit of a radiation-amplitude vacuum-like if and exactly if the quark line 
of the leading final state quark (or one of the two quark and anti-quark lines of the most energetic final state gluon) 
connect to the radiated gluon without passing through components of the medium. 
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then provide probability-conserving corrections to the same color configuration. Such corrections 
contain virtual terms, as exemplified by the N = 2 opacity contribution to [aai]. 

In general, to each order in opacity, only one vacuum-like color configuration [ 
opens up anew. However, with increasing order in opacity, more and more medium-modified color 
configurations can arise. For instance, up to second order in opacity, we find the terms 

[ai a] , [a2 a ai] , [a2 ai a] , [ai a 02] , [ai 02 0] ■ (4.27) 

The strong increase in the number of medium-modified color configurations with increasing order 
in opacity, and the corresponding decreasing weight of vacuum-like contributions (as seen in the 
explicit N = 2 calculation) both support the physically intuitive idea that with increasing density 
of the medium it becomes easier for a radiated gluon to decohere in color from its partonic sister 
fragments. Our calculations show that the contribution of color-decohered gluons to the medium- 
induced radiation always exceeds 50 % whenever the latter is sizable, and it may well be much 
larger than 50 %. 

5 Hadronizing parton showers with medium-modified color fiow 

Since QCD conserves color, a dynamically consistent hadronization model must respect color fiow in 
interfacing the fragments of perturbative parton showers with the hadronic final states. Thus, even if 
hadronization occurs time-delayed and outside the medium, it can be affected by the modification 
of the color connections introduced by the medium. The Lund string-fragmentation model is a 
prime example of a phenomenologically successful hadronization model that respects color fiow. 
It is very well documented in the literature [18, 21]. Hadronization is modeled by decomposing 
each parton shower into a set of color singlet 'strings'. These strings stretch between q and q 
endpoints and represent color-connected gluons as kinks. They are then decayed into hadrons 
through the excitation of qq pairs from the vacuum. In section 2 we already discussed some examples 
of interfacing the final stage of a branching process with LuND strings. Since the medium-modified 
parton branchings studied in sections 3 and 4 define color fiow unambiguously, the above string- 
fragmentation routine can be applied without additional model assumptions. In this section, we 
explore how the distribution of hadronic fragments can be affected by medium-modifications of color 
fiow. To this end, we interface the possible color-configurations arising from the medium-induced 
branching of high-energy partons in the plasma with the LuND hadronization model. Numerical 
results were obtained with the string-fragmentation routine implemented in Pythia 6.4 [18]. The 
strings to be decayed were built using the routines PYIENT (to add a given parton to the event) 
and FY JOIN (to join the partons provided by the user into a string, according to the proper order). 
The aim of the following analysis is to illustrate, for specific examples, how the distribution of 
hadronic fragments can be affected by a medium-modified color fiow. 

We start by considering Fig. 12. To first order in opacity, medium-induced gluon radiation of a 
quark projectile exhibits two different color flows. In the channel labeled as Final State Radiation, 
the emitted gluon will be part of the same Lund string as the leading quark. In contrast, in the 
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Figure 12. The = 1 opacity correction to the fragmentation of a liigh-pT quark, interfaced with the 
LuND-string hadronization routine. In the case of final state radiation (FSR, left hand side) the radiated 
gluon remains color connected with the other daughter of the branching and, if sufficiently collinear, can 
contribute to the energy of the leading hadron. For initial state radiation (ISR, right hand side), the 
gluon results color decohered from the leading projectile fragment, independently on the emission angle; 
an independent Lund string is associated to it, whose decay will contribute to an enhanced multiplicity of 
soft particles. 

case of Initial State Radiation (right panel of Fig. 12), the radiated gluon is color decohered from 
the projectile fragment; this means that it is instead linked through an independent string to a 
low-pT particle in the medium. For the following, the Lund strings in Fig. 12 are defined in terms 
of the 4-momenta of their end-points and kinks. For simplicity both the leading quark fragment and 
the radiated gluon are assumed to be emitted at mid-rapidity (r/ = 0, i.e. 6 = 7r/2) and at relative 
azimuthal angle 0. The other endpoint of the string will then be attached either to a particle from 
the medium or from the beam remnant: in both cases it will sit at low-p-r- Medium particles are 
taken with a typical thermal energy ~ 3T, with random event-by-event momentum orientation. For 
the temperature we take the value T = 200 MeV^. For the antiquark from the beam remnant, we 
choose a large momentum along the beam direction {E = 1 TeV and pt = 0)- The only dependence 
on E resides in the extension of the plateau in the rapidity distribution of (very soft) hadrons from 
the decay of the subleading string. 

Figure 13 shows examples of the distribution of hadronic fragments of the Lund strings de- 
picted in Fig. 12. These results were obtained for a typical partonic configuration with a quark 
at high transverse momentum pj- and a radiated gluon at much smaller transverse momentum (as 
an illustration we chose here and in the following fc^ = 0.1 pt)- In the left panel of Fig. 13 we 
display, for the two different color channels, the fragmentation pattern of a hard quark branching 
in the medium. The FSR case (red curve) corresponds to a vacuum-like color flow: in this case 

^The results were found to have a negligible dependence on tlie temperature. This is consistent with the idea that 
the precise position of the soft endpoint of a string is unimportant for the hadronization of hard fragments. 
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Figure 13. The px and rj distributions of the hadrons from the fragmentation of the Lund strings shown 
in Fig. 12. Both the quark and the gluon are emitted at midrapidity at relative angle (p = 0.1. Left 
panel: fragmentation pattern in the FSR (in red) and ISR (in green) color channels. Right panel: rapidity 
distribution of the hadrons in the ISR channel. The sharpest peak around to r/ = (continuous line) comes 
from the fragmentation of the leading string. The pattern "broad peak + plateau" (dashed line) arises 
from the fragmentation of the subleading string, connected to the beam remnant (hence the long plateau). 
Also shown (dot-dashed line) is the case in which both endpoints of the subleading string are attached to 
a medium particle. 



there is hadronic yield in a transverse momentum range that exceeds the px of the leading quark. 
In the Lund model, this accounts for the fact that QCD is a finite resolution theory in which a 
perturbatively radiated gluon does not automatically increase the hadronic multiplicity by order 
unity or more: it is not necessarily 'lost' but, remaining color-connected with the other daughter of 
the branching, may still contribute to the formation of the leading hadron. In contrast, the ISR case 
(green curve) clearly shows that medium modification of color connections between the radiated 
gluon and the projectile fragment results in a softening of the hadron distribution: all hadronic 
yield above px is suppressed and an additional contribution arises at soft momenta below kx- The 
reason is that, for the ISR contribution, the color- decohered gluon and quark belong to different 
strings and thus cannot contribute to the same leading hadronic fragment. Therefore, hadronic 
multiplicity increases by construction with each color- decohered gluon by order unity or more, and 
the additional multiplicity is found in soft fragments of transverse momentum lower than kx, which 
is much smaller than px- 

These differences in the color flow of the ISR and FSR contribution have consequences for 
the distribution of hadronic fragments. In particular, the fragmentation of the LUND string of 
a vacuum-like (FSR) contribution results mainly in semi-hard and hard hadrons. For instance, 
fragmentation of the FSR string of total energy ~ 55 GeV in Fig. 13 yields on average (Nh) = 5.4 
hadrons, of which 3.9 carry px > 2 GeV transverse momentum. Since the multiplicity of Lund 
strings grows only mildly with the total length and with the number of small kinks, the string 
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of the ISR contribution that contains the leading quark fragment will decay into almost as many 
hadrons {{Nh) = 5.2 for the case shown in in Fig. 13). However, these hadrons carry a smaller total 
transverse momentum, the remaining fraction kx / {kj- + pj-) being carried by the fragments of the 
subleading string. The latter stretches over a short distance in transverse momentum, but it can 
stretch over a long distance in rapidity, and hence it can yield high multiplicity. For instance, for 
one of the cases shown in Fig. 13, the string stretches from a medium component around t] = 
to a beam remnant at projectile rapidity and thus distributes (^h)^bi(Jad~''°^™ ~ hadrons over 
approximately 10 units in rapidity. Alternatively, it is conceivable that due to further interactions 
between the medium and the radiated gluon, both ends of the subleading string connect to soft 
components close to mid-rapidity. In this case, the subleading string produces a much smaller 
number of hadrons {{Nh)'^^i~^^'^'^'^ = 3.7 for the case in Fig. 13); the hadronic distribution close to 
the rapidity ?7 = of the jet, however, changes only mildly. 

Recent measurements by the CMS collaboration indicate that the medium-modified fragmen- 
tation of jets in heavy ion collisions at the LHC is accompanied by a significant increase in soft 
hadronic multiplicity ^ that appears outside typical jet cones at large angle with respect to the jet 
axis [4], while the hard part of the jet fragmentation function is consistent with that of a vacuum 
jet of lower transverse energy [24]. Here, we note that medium- modified color decoherence is a 
natural candidate mechanism for two of these qualitative features. Namely, the leading string of 
the ISR contribution coincides by construction with that of a vacuum quark jet of lower transverse 
energy {pt, rather than pt + ^t)- Moreover, medium- induced color decoherence leads naturally 
to distributing the transverse energy kx amongst very soft components. It is less clear, however, 
whether this additional soft multiplicity is distributed naturally over large angles outside the jet 
cone. Our study indicates that the color decoherence associated to the = 1 opacity calculation 
is not enough to achieve such an angular broadening and a larger amount of rescattering at the 
partonic level (that can be only achieved going to higher orders in opacity) is mandatory. We cau- 
tion that the numerical studies presented here focused entirely on the effects of color decoherence 
while a phenomenologically satisfactory description of jet quenching requires of course to interface 
these effects with the medium- modification of kinematic distributions. Consistent with previous 
discussions we therefore expect that the broadening of soft components is mainly due to kinematic 
effects (see Refs. [25-28] for models consistent with this idea). 

We finally comment on the potential relevance of our findings for single-hadron spectra which are 
particularly sensitive to the modifications occurring to the hard tail of the fragmentation function. 
It is customary to characterize this effect by considering a steeply-falling parton spectrum 1/p^ 
that if convoluted with a Fragmentation Function (FF) yields dcr^^'^'^°^ ~ {x"'~^)/p^. This shows 
that single hadron spectra are mainly sensitive to a higher {{n — 1)*^) moment of the FF. In 
Fig. 13, we display the ratio of the 6-th moments for different color channels in a specific kinematic 
configuration of a parton shower. This figure illustrates that for an identical kinematic partonic 
configuration (i.e., for the same amount of energy loss at the partonic level), the modification of 

^An increase in soft hadronic multiplicity is also seen in hadron- and photon-triggered jet-like distributions at 
RHIC. For recent preliminary data, see e.g. Refs. [22. 23]. 
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Figure 14. The ratio of 6-th moments of the ISR and FSR contributions to the fragmentation functions 
of partonic configurations shown in Fig. 12. Here, prior to hadronization the leading quark has transverse 
momentum pT, and the radiated gluon carries a scaled transverse momentum kg = 0.1 pr and is emitted 
at (j) = 0.1. 



the color connections can introduce an additional significant source of suppression of the hadronic 
spectrum. Quantitatively, in the case displayed in Fig. 13, one finds for instance (x^)isr = 0.052 and 
{x^)fsr = 0.078 for the ISR and FSR channels respectively. In summary: since medium- induced 
color decoherence of radiated gluons depletes naturally the high-p^p tail of fragmentation functions, 
it can contribute to reducing the nuclear modification factor -Raa- A similar observation was made 
in a recent study [19], in which parton splittings with medium-modified color flow were interfaced 
with a cluster hadronization model. Also in that case the suppression showed a weak dependence 
on transverse momenta and persisted at transverse momenta exceeding 100 GeV. Fig. 14 supports 
the picture that medium- induced color decoherence can be a relevant factor for understanding the 
value and pT-dependence of the nuclear modification factor at high transverse momentum. 

Elementary formation time arguments suggest [29] that hadronization can occur time-delayed 
outside the medium, and that the earliest parton branchings in a shower can be medium-modified 
since they occur inside the medium. However, these earliest branchings may be followed by subse- 
quent partonic branchings outside the medium prior to hadronization. The question arises whether 
the features observed in Fig. 13 and 14 are robust against inclusion of such subsequent parton 
branchings. The structure of typical parton branching histories indicates that inclusion of a single 
additional splitting on top of the ISR and FSR contributions of Fig. 12 provides already significant 
insight into the question of how final state branchings may alter the effects of medium- induced color 
decoherence on hadronization ^. As shown in Fig. 16, inclusion of another parton splitting adds an 

^ As a consequence of the steeply falling px dependence of spectra, by triggering on the px of a parton at the end 
of a shower evolution one typically selects parton showers that have suffered a very limited number of branching. 
For instance, for a sample of quarks at LHC energy generated by Pythia in a 'hard event' + 'final state shower' 
process, one finds for a trigger px ^ 50 GeV of the final quark an average number of daughter gluons of order ^ 2.5. 
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Figure 15. The same as Fig. 13, but for two choices of quark energy, and for a varying azimuthal angle 
(j) of the radiated gluon. As the relative angle gets larger the FF suffers a softening. 




Figure 16. The color flow associated with medium-induced FSR (left hand side) and ISR (right hand 
side) in the = 1 opacity expansion, supplemented with an additional second gluon emission at late time. 



additional kink to the leading string, but it does not change the fact that the color- decohered gluon 
of the ISR configuration hadronizes in a separate, subleading string, while in the FSR contribution 
both gluons are part of the leading string. As a consequence, if both gluons of the FSR contribution 
are produced sufficiently collinear, they can contribute to the production of the leading hadron, and 
the corresponding fragmentation function in Fig. 17 extends up to the transverse momenta that 
correspond to the sum of the transverse momentum of the three projectile partons in the final 

This number increases only very weakly with p^'^. 
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state. In contrast, the color- decohered gluon of the ISR contribution is lost for the formation of 
the leading hadron. In this way, Figure 17 clearly illustrates that additional parton branching at 
late times does not wash out the effects of medium-induced color decoherence. More precisely, 
medium- induced color decoherence removes efficiently a fraction of a jet's energy from the high-pr 
part of the fragmentation function by hadronizing it independently in soft fragments. The leading 
string, however, irrespective of the number of gluons (i.e. kinks) radiated outside the medium that 
it involves, appears to correspond at sufficiently high px to a vacuum fragmentation pattern of a jet 
of lower (ISR gluons are decohered) total transverse energy. We emphasize that this observation 
depends mainly on two generic features, namely that color-decohered gluons correspond to different 
strings and that additional branchings at late times correspond to (typically small) kinks on the 
same string and thus contribute to the distribution of leading fragments. 




p.j. primary hadrons (GeV) 



Figure 17. The pr-distribution of hadrons from the fragmentation of the Lund strings interfaced with 
the two color-flows involved at = 1 order in opacity supplemented by the subsequent showering shown 
in Fig. 16. The two parton branchings were at azimuthal angle (j) = 0-2 and (j) = 0.1, respectively. The 
figure illustrates that the effects of medium-induced color decoherence are not washed out by subsequent 
showering. 



6 Conclusions 

In the absence of medium effects, the only color correlation between a final state parton shower 
and the rest of a hadronic collision is the one between the color of the primary 'parent' parton and 
some color- compensating beam remnant. Jet-medium interactions change this feature generically 
and characteristically. To date, very few studies [14, 15, 30-35] have addressed aspects of the role 
played by medium modifications of color flow. Most studies focus solely on medium- induced changes 
of the multiplicity and kinematic distribution of partons in the shower. In the present work, we 
have pointed out that jet-medium interactions give rise to qualitatively novel color correlations in 
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which components of the parton shower are color- decohered from their sister partons in branching 
processes. 

Let us summarize our findings in a formulation that is slightly different from the one adopted so 
far. Physically, hadronization of parton showers amounts to a procedure that specifies the overlap 
between the color singlet components in the shower and the physical Hilbert space of all hadronic 
wave functions corresponding to experimentally accessible states. As is generally known and as we 
have illustrated for the Lund hadronization model in Fig. 15, this mapping can result in hadronic 
fragments that are more energetic than the most energetic partonic component that is hadronized. 
This is so, since a radiated partonic component - as long as not decohered completely from its sisters 
in the branching process - can display an overlap with the same hadronic one-particle state as its 
sisters. In this sense, the energy radiated perturbatively in gluons is not automatically an energy 
lost for the formation of a hadron; it is only 'lost' if the gluon is decohered. In the vacuum, sister 
partons typically remain color connected and their decoherence therefore arises only kinematically, 
as a consequence of their relative transverse separation. To date, studies of parton energy loss 
have mainly focused on modifications of this kinematic decoherence. Here, we have demonstrated 
that in comparison to purely kinematic effects, medium-modified color flow is a highly efficient 
decoherence mechanism and can thus contribute significantly to softening hadronic distributions 
(see section 5). As emphasized in section 2 this is a qualitatively novel and generic property of 
parton energy loss mechanisms that arises from the gluon exchanges between the parton shower 
and the medium, independently of model-specific dynamical details. It will generically increase the 
number of decohered color singlets in the parton shower. The relative rate of partonic configurations 
with an increased number of color-decohered components is large: in the explicit color-differential 
calculations presented in sections 3 and 4, it always exceeds 50 % of the entire medium-modified 
gluon radiation. 

Remarkably, we found in section 5 that the medium-induced color decoherence of gluons from 
the parton shower leads naturally to their hadronization into soft components while the most ener- 
getic components will hadronize like vacuum structures of reduced transverse energy. This appears 
to be characteristically different from the typically studied kinematically induced medium modifi- 
cations where one expects generically that the enhancement of soft jet fragments is accompanied 
by medium-modifications of the shape of the fragmentation function at all momentum scales. On 
the other hand, it is qualitatively in line with CMS data on dijet asymmetries. In particular, CMS 
data show jet fragmentation functions that are at high transverse momentum consistent with the 
fragmentation of vacuum jets of degraded total transverse energy Et — Ai?. Moreover, the miss- 
ing jet energy AE is recovered in very soft hadrons. Both these features arise naturally from the 
medium-modified color decoherence of jet fragments studied here. We caution, however, that the 
experimentally observed soft components are distributed over a wide region in Arj x A0. While we 
have seen some broadening of the distribution of soft fragments (see Fig. 15), we did not identify a 
generic argument that could account for such an observed wide distributions in the presence of only 
= 1 gluon exchange with the medium. It remains to be clarified for instance to what extent soft 
fragments are broadened further if subsequent medium-induced gluon branchings or higher orders 
in opacity are taken into account. More generally, the present work was limited to showing that 
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medium-induced color decoherence is a non-negotiable aspect of parton energy loss mechanisms 
that can lead to qualitatively novel and numerically significant features in jet quenching calcula- 
tions. The next step is now to incorporate the generic features of medium-induced color decoherence 
established here in a more complete dynamical modeling of jet quenching. 

Acknowledgements 

We thank N. Armesto, C. Salgado, B.G. Zakharov and K. Zapp for useful discussions. JGM 
acknowledges the support of Fundagao para a Ciencia e a Tecnologia (Portugal) under project 
CERN/FP/116379/2010. 

A Gluon projectile: color-differential spectrum 

Here we provide further details on the color-differential = 1 opacity calculation for the case of 
a projectile gluon impacting from the distant past. This extends the discussion of the diagrams of 
Fig. 6 in section 3.1 to the case of a medium- modified splitting 



As indicated above, the incoming gluon has color d. The outgoing gluons carry momentum fractions 
Xg and ^ — Xg with Xg <^ 1 — Xg. The color b is attributed to the outgoing harder gluon that carries 
momentum fraction 1 — Xg and to which we refer as outgoing 'projectile' gluon; the softer outgoing 
gluon carries color a and the gluon exchanged with the medium has color index ai. The polarization 
vectors for the incoming and outgoing 'projectile' gluon, as well as for the radiated gluon are 



respectively. For a high energy gluon, the interaction with the medium is described by the eikonal 
vertex gf°''"^g'^P{p + p'y. The amplitude A^(a) of a gluon that first interacts with the medium prior 
to gluon radiation reads then 
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In close analogy, we find 
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and 

^ -1 
^M^c) = Y.^-i)9'n'f'\e,.ef)——-^2{l-x,) [e,- (fc,-q)]2p+^(q)e^''-T- . (A.5) 

n=l ^'^S Q) 

These expressions are given in terms of the generators of the adjoint representation, (Tl)^,^ = if"'^'^- 
For a discussion of color flow in the large- iVc limit it is useful to express them in terms of traces 
over products of generators in the fundamental representation. For instance, starting from the 
normahzation Tr(t"t^) = 1 /25"^ and if'''"' = 2 Tr ( [f, t''] t") , one finds 

(T'*T"i)m = 2Tr ([^^r][r^t'^]) . (A.6) 

As explained in section 3.2, this allows one to separate the entire gluon radiation amplitude into 
six different color channels. The corresponding diagrammatic contributions are shown in Fig 5 
and listed in (3.15). Here, we provide their explicit expressions. Representing the incoming gluon 
projectile in the large- A'^c limit by a quark and an anti-quark leg, we label these contributions as 
follows: 

• FSR(g): final state radiation of the quark leg (configuration haaid) 

^M,..,, = j:{-^)2Tr[thH^H^) [ "[l^'^r.'jy ' ^ff^) 2(l-x,)(e..,).M-'^"^ (A.7) 

• FSR(g): final state radiation of the anti-quark leg (configuration hdaio) 

rMua.a = B-.) 2Tr (nVt-) (^^^^^ - 7^^) 2(l-x,)(e..,)Al-'(") (A.8) 

• ISR(g): initial state radiation of the quark leg (configuration haiad) 



N 



ISR(g): initial state radiation of the anti-quark leg (configuration bdaai 



^-M..... = E -(-0 2Tr {tH'^HH^) f ^ - ^^^^ ) 2(l-x,)(6,.6,)>l-'(") (A.9) 



^A^...a. = g -H) 2Tr (t-tV^) - ^^^^ j 2(l-x,)(6,.6,)>l-'(") (A.IO) 

RqSg: radiation of the quark leg but medium interaction on anti-quark leg, configuration 

(badai) 

iA1.a.a.=f:-(-02Tr(tWt-) f - ^) 2(l-^.)(^re/).M:/'^") (A.ll) 

[Kg XgQ) fS 
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• RgSq: radiation of the anti-quark leg but medium interaction on quark leg (configuration 
baida) 

^^^..... = f:-(-02Tr(rt-^^) - ^] 2{l-x,){e,e,)MT^-^ (A.12) 

In the above, the first four channels are consistent with the simplified large- A^^^ picture, where the 
radiation from a gluon is viewed as the incoherent superposition of the medium-induced radiation 
from its quark and anti-quark legs. The last two channels, on the other hand, represent contributions 
for which the q (q) leg radiates and the q (g) leg exchanges color with the medium. For these 
contributions, a distinction between initial state and final state radiation is not needed, since it 
does not affect the color ordering of the amplitude. However, in the soft ^ 1 limit, these last 
two channels give a subleading contribution. Squaring these two amplitudes one would find a QED- 
like ('Bethe-Heitler') radiation spectrum that is 0{x^). On the other hand, the first four channels 
are not suppressed in the soft limit and provide, to leading order in Xg, an identical contribution: 
they sum up to the well-known Gunion-Bertsch spectrum. 



B Details on the N = 2 calculation 

In section 4 we quoted the result for the color-differential radiation spectrum at order = 2 in 
opacity. The physical meaning of the various contributions became particularly transparent by 
expressing the result in terms of the effective currents Jj (z = 1, 2, 12). Here we give details on how 
to derive these results. For this purpose, the calculation of a few diagrams will be sufficient, the 
other ones will follow from the modified dependence on the currents Jj. 

B.l Direct contribution 



Figure 18. A contribution to with its corresponding '02010' color flow (which gets contribution also 
from the diagram with the radiated gluon scattering on the first scattering center). 



We start by considering the direct term in the N = 2 spectrum. For illustration, we focus on 
the 'a2aia' channel that can be written as a sum of two terms, 

Mr^'^^Mgl^ + tMii'l,, (B.l) 
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M[a_2a_la] ' \^ M*[a_2a_la] 

Figure 19. The color flow in {\Ai2^"'^"'\'^). Each disconnected loop gives rise to a factor Nc- Each gluon 
line provides a factor 1/2 (due to the Fierz identity). 



where A^ijaia is the contribution with the gluon radiated before any interaction with the medium 
see Fig. 18) and Mi['Xa describes the radiated gluon interacting with the first scattering center. 
The first term takes the explicit form 

'■^'^a.a.a- 2^9^^ ^ M I ^ („,) ^ („,) ^ J 271 J {2n)^ J {27r)^ 

X e^(«r+92"Xe^''2-«-0(^)[(2p^ - q2)+A{q2)]- 



{pf-q2)^-M^ + iT] 



X m2pj - 2q, - gi)+^fa)] ^^^_^^^, '_^,_^^^ (0[(2p. - k,) . e,]. (B.2) 



Here, the pole of the propagator {Q = qi + q2) 



{Pf - q2? - + zr/ 2(1 - x,)p+[q^ - (^^^^ffeg^ + iv)] 
allows one to perform the integration. Furthermore, momentum conservation implies 
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Up to 0{xg) corrections, one gets 
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X e-°-n A{q,) Aiq^) | 2e, ■ ^^^^ ] 2p+. (B.5) 
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Proceeding analogously with the contribution M.a2aia-, one obtains 



{2'nYJ (27r)2 



ni,n2 

X e*^°""i^(gi)^(q2) (2e, ■ Ji) 2p+. (B.6) 

After squaring, tracing over colors 

— Tr ( T"-' T"'' ^ — —() ,T r^"'"'' CR T^l 
and averaging the phase factor over the transverse location of the scattering centers 

I rf^„e-('^"-'^"')-^ = -^(2vr)25(g„ - , (B.8) 

one finds 



which can be written as 

/ ^cl \ 2 / -^ ^^el \ /• /I ^^cl \ 

|Ji|2(2p+)2. (B.IO) 



ni,7i2 ^ ' - ^ -J- ^ . \ 

Here, cx''^ denotes the quark elastic cross section. One has still to perform an average over the 

longitudinal position of the scattering centers. Notice that in the above we assumed x^^ > x^_^ 
which implies 



L+ J J. "1 L+ V L^ 

ni,n2 ^ / >' u >' 



Recasting this expression in the form 



one finds 



Ai J 2 \AiL+ J 2 VA+/2 



(B.12) 



which coincides with the expression quoted in the text. Finally, the longitudinal average of terms 
containing interference factors (which appear in other channels) can be performed exploiting 

E (^) ' r ^< /r ^< - = (^-^^^ 

ni,n2 ^ / -^0 Jxt-^ { ) 
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The medium average Eq (B.7) is the same as employed in standard color inclusive calculations. 
It traces independently over the color of the scattering centers ni and n2. We emphasize that in 
the large- Ai'c limit this procedure does not affect the decomposition of the radiation spectrum into 
distinct color channels. Here, we consider for simplicity the case in which the scattering of the 
projectile occurs on a medium particle in the ( ant i-) fundamental representation. After averaging 
over all the incoming colors, one has, according to Eq. (B.9), the overall color factor 

^^^^ - ^^2iV^2iV^ ^ 32 • ^^-^^^ 

This result can also be obtained graphically from Fig. 19 with the simple rule that each closed loop 
produces a factor A'c, and each gluon exchange, due to the Fierz identity 



^Ij^i = \^ii^jk - -^^ijhi ~ ]j^5ii5jk , (B.16) 



gives rise in the large- A'^c limit to a factor 1/2. One gets then for Fig. 19, after averaging over the 
initial colors: 

in agreement with Eq. (B.15). 
B.2 Virtual contribution 

Let us now address the evaluation of the 2Re(A^™*A^*) contribution. The amplitude Ai^ interferes 
either with the initial-state radiation (oja channel) 

zMT = ^Mr = {^g) Y.it'^H'') y^e-^'^'-^^ 2 e,- J(A;„ q') e^^"^^ (2p+) A{q') , (B.18) 

or with the final-state radiation (aai channel) 

iMl'' ^ iMT = -iw) Y^^tH^') /(^^"""'"^ 2 e,-J{k,, q') e^^°^^ (2p+) ^(q') Tg . (B.19) 

The various terms contributing to A^g'''* can be obtained from the A^ = 3 diagrams listed in [36]. 
While each individual diagram can produce up to four different contributions (xi =X2 or X2 =2^3, 
interference with /A^^ and A4.\^), only some of them yield a leading contribution in the large- Ac 
limit. Their expression can be directly read from [36] after taking properly the contact limit of two 
of the three scattering centers. To fix the exact overall factors it is sufficient to consider one case, 
e.g. the term with the gluon radiated before the interaction with the first scattering center and 
with a double interaction at the second one (see Fig. 11). This corresponds to the first of the 15 
diagrams shown in Appendix C of Ref. [36] (hence we label it as '/'). After setting X2 = x^ and 
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considering the projectile arriving on- shell, one has 

dqi f dq^ f dq^ [dq^^ fdq-^ 

i<j 



(27r)2 J (27r)2 J (27r)2 j 2ii j 2^ 



X e 



m2pf-q,)+A{q,)] 



i)[{2pf-2q^-q2yA{q^)] 



{Pf - 93)^ - M'^ + if] 

x{i)[{2pJ-2qs-2q2-qlyA{q^)]- 



{Pf -q2- qzY - M'^ + ir] 



m2p,~k,)-e,], (B.20) 



(pi - k^y - + tr]' 

where ^23 = '?2^+5'3'- Furthermore, from the on-shell condition and momentum conservation one has 



(?i+?2+g3) 



2XgP+ 2{1-Xg)p+ 2p^ 



{Q = Qi + q2 + Q3) 



(B.21] 



This fixes q-^ . The integration over can be evaluated picking the pole in the complex upper 
half-plane of the propagator 



{pf - q2 - qsY - + i7] 2{1-Xg)p+ 



^23 



(B.22) 



which in the high-energy limit lies at ~ it]. For what concerns the integration over q^ , in the 
contact limit, no phase factor is present and one no longer closes the contour in the complex plane. 
However one can rely on the following representation of the Dirac delta (no contribution to the 
integration arises from the imaginary part): 



{pf - qsY - + tr] 2{1-Xg)p^ 



IT 



2(l-x,)p+ + *V 



2{1-Xg)p+ TT 

We thus obtain in the soft limit 
/ ■^9 



^3 2il-Xg)p+ 



2(1 - Xg)p + 



5iq^). (B.23) 



• « /(virt (I) 



i<j 



dqi f dq^ f dq: 



(27r)2 J (27r)2 J {2tt) 



X e-''^i "'e-*('^2+'^3)-"^e*"o"«^^(g3)^(q2) {2p+) A^q^) \2eg-^ 



k,. 



(B.24) 



We now evaluate the interference with the = 1 amplitude, taking as usual the proper medium av- 
erage (over colors, transverse and longitudinal location of the scattering centers). The contribution 
from the final-state radiation channel is vanishing. 



2Re(Ai:;X2.(K'T) = 0, 



(B.25) 
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For the interference with the initial state radiation one gets 



2R<=(K";::jA<n->=-49^c.(^;)^j_.i^i<J<J^^^ . (b.26) 

The other contributions can be evaluated analogously. 
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